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ABSTRACT 

The solar surface is covered by high-speed jets transporting mass and energy 
into the solar corona and feeding the solar wind. The most prominent of these 
jets have been known as spicules. However, the mechanism initiating these erup- 
tions events is still unknown. Using realistic numerical simulations we find that 
small-scale eruptions are produced by ubiquitous magnetized vortex tubes gener- 
ated by the Sun's turbulent convection in subsurface layers. The swirling vortex 
tubes (resembling tornadoes) penetrate into the solar atmosphere, capture and 
stretch background magnetic field, and push surrounding material up, generating 
quasiperiodic shocks. Our simulations reveal a complicated high-speed flow pat- 
terns, and thermodynamic and magnetic structure in the erupting vortex tubes. 
We found that the eruptions are initiated in the subsurface layers and are driven 
by the high-pressure gradients in the subphotosphere and photosphere, and by 
the Lorentz force in the higher atmosphere layers. 

Subject headings: Sun: photosphere, chromosphere, magnetic fields; Methods: 
numerical; MHD, plasmas, turbulence 

1. Introduction 

One of the most frequently observed phenomena of dynamical coupling between the 
solar convection and atmospheric layers are plasma eruptions on different scales, such as 
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CME, flares and spicules. Rapid increa se of the observ ational power with new instruments 



e.g. HMI/SDO flScherrer et al. I l2012h. SOT/Hinode flTsuneta et al. I l2008h. NST/BBSO 



jGoode et al~ll2010h . IMAX/SUNRISE JSolanki et al. I boioh allows us to detect various 



features on smaller and smaller scales, ans measure their properties. Realistic numerical sim- 



ulations, based on 'in intio 1 phy sical principles (e.g lStein fc Nordlund Ill998l ; IJacoutot et al. 
2008at iKitiashvili et al. Il2012al ). have been able to reproduce and understand physics of the 
observed phenomena, and predict new effects, which are hard to detect in observations. In 
this paper, we present new results of radiative MHD simulations, which shed light on the 
mechanism of small-scale eruptions in the solar atmosphere, and link it to the dynamics of 
turbulent magnetized vortex tubes. 

Previously it was shown that an important role in dynamics the turbulent convective 
layers and low atmosp here, and in the chromo s phere heating is played by small - scale vortex 



tubes , by simulations ( iKitiashvili et al. 112010 



2012 



Balmaceda et al. 2010; Cap et al. 2010; Yurchyshyn et al. 



ispnerj 
. 2011 



2012a; 



Steiner et al. 



2010; 



Moll et al. 



2011 ) and high-resolution observations (e . g. IWedemever-Bohm fc Rouppe van der Voort 112009 



2011 



Fedun et al. 112011; Ji et al. 



Wedemeyer-Bohm et al. 1 120121 ) . In particular, the numerical simulations revealed the 



dynamics of the vortex tubes; and it was found that despite preferentially strong downflows 
in the vortex core, strong upflows may occur quasiper i odicall y, pushing matherial into the 
atmospheric layers (see in Fig. 5 in IKitiashvili et al. I ( I2012a[ )). In this paper, we demon- 
strate that these quasi-periodic upflows associated with the vortex tube dynamics, in the 
presence of magnetic field produce small-scale jet-like ejections of plasma generating shocks 
in the atmosphere, and discuss the origin and properties of these eruptions. 



2. Computational setup 



To simulate the coupled dynamics of top layers of the turbulent convective zone, pho- 
tosphere and low atmosphere of the Sun we use the 3D radiative MHD 'Solar Box' code 
deve loped at the Stanford C enter for Turbulence Research and NASA Ames Research Cen- 



ter ( IJacoutot et al. Il2008al ). The code impl ements LES turbule nce models, the real-gas 



equation of state, astrophysical opacity tabl es ([Rogers et al. Ill996l). and the sta ndard model 



of the solar interior for initial conditions (jChristensen-Dalsgaard et al. II1996I ). Radiative 



transfer between fluid elements is calculated using a 3D multi-spectral-bin method, and a 
long-characteristic approach, assuming the local thermodynamic equilibrium. 

The physical description of dynamical properties of the solar convection is improved 
through the implementation of subgrid-scale LES turbulence models, which can effectively 
increase the Reynolds number and provide representation of small-scale motions closer to 
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form flMoin et al. Nl991[ ): r, 



the reality (IJacoutot et al. Il2008al ). Here we used a Smagorinsky eddy- viscosity model 



(ISmagorinskv 1 119631). in which the compressible Reynolds stresses were calculated in the 



-2C s A 2 \S\(S id - u^kSij/3) + 2C c A 2 \S\ 2 5 ij /3, where the 



Smagorinsky coefficients are C$ = Cq = 0.001, Sij is the large-scale stress tensor, 
A = (dx x dy x dz) 1 ' 3 with dx, dy and dz being the grid-cell dimensions. 



and 



The presented simulation results are obtained for a computational domain: 6.4 x 
6.4 x 5.5 Mm 3 (with a 1-Mm layer above the photosphere, and a 5-Mm layer below) with 
12. 5 2 x 10 km grid-cells. We modeling the conditions of a quiet-Sun region: with an initially 
introduced vertical magnetic field, B z q = 10 G, and, for the comparison for pure hydro- 
dynamic case. The lateral boundary conditions are periodic. The top boundary is open to 
mass, momentum and energy transfers, and also to the radiation flux. The bottom boundary 
is open only for radiation, and simulates the energy input from the inter ior of the Sun. The 



simula tion results previously were compared with a similar type c ode by iNordlund fc Stein 



(200 11) for some test ca ses, and with photospheric observations ( IJacoutot et al. I l2008aU b 



Kitiashvili et al. 1 120131 ) . Extending the computational box domain into the atmosphere 



allows us to model effects of the intensive energy exchange between the photosphere and 
the chromosphere, and investigate some of most the energetic phenomena in the quiet Sun: 
spontaneous high-speed flow eruptions along the magnetic flux tubes. 



3. Flow ejection due to vortex tube dynamics 

The highly-turbulent solar subsurface layers are a place of the origin of numerous vortex 
tubes, formed due to the granular overturning and Kelvin-Helmholtz instabilities (see Kiti- 
ashvili et al. 2012 and references therein). The turbulent vortex tubes are usually located 
in the intergranular lanes; they often become vertically oriented, penetrate above the solar 
surface, and can be stable for longer than a typical granulation life-time. High-speed, near 
sonic, flows of the vortex tubes are accompanying by sharp variations of temperature, density 
and gas pressure, and can strongly affect the dynamics of their environment. According to 
the previous numerical simulations, the vortex tubes penetrate into the chromospheric layers 
in both, magnetic and nonmagnetic cases. In the quiet-Sun region, the magnetic field effects 
reveal themselves mostly in the higher atmosphere as magnification of the hydrodynamic 
effects playing the dominant role in the turbulent surface and subsurface layers. Thus, the 
contribution of the magnetized vortex tubes to the chromospheric dynamic s and energetics 



is more significant than in the hydrodynamic case (jKitiashvili et al. Il2012al ) 



Our simulations show that the upward vortex tube penetration into the solar atmosphere 
is often quasi-periodic and accompanied by spontaneous flow ejections. A time-sequence of 



-4- 



various vortex tube properties with a cadence of 15 sec for a strong event is illustrated in 
Figure [U where panels a — c) show horizontal snapshots of temperature, vertical velocity 
and density at height h = 625 km; and panel d) shows vertical cuts of log(p) taken along 
the x-axis through the region marked by short white lines in panels a — c) . Figure [T] shows a 
complicated structure and dynamics of the eruptions with downflows in the vortex core, and 
the upward eruption flows in the surrounding region. As seen in the temperature distribution 
plots in panel a) these eruptions are hotter than the surrounding plasma, and can provide an 
extra chromospheric heating in addition to the heating through the vortex core. The shape 
of the flow ejections in the vortex tubes is not completely circular due to the non-vertical 
vortex tube orientation, and interaction with surrounding flows (Fig. [1] c, d). 

The vortex tube causes strong swirling motions of the subsurface and atmospheric layers 
around the vortex core. Figure [2] illustrates various stages of the flow ejection, where the 
streamlines show a general behavior of the velocity field; yellow-blue isosurfaces represent 
the pressure gradient, normalized by density, of 5 x 10 4 cm/s 2 (yellow) and —5 x 10 4 cm/s 2 
(blue); and the grey isosurface shows the temperature distribution of T = 6400 K. In the 
atmospheric layers, the flow electio n starts with the process of formation of a vertically 



oriented vortex tube, as described by lKitiashvili et al. I (J2012bJ). This creates strong vertical 



pressure gradients: negative in the vortex core and positive at the vortex periphery (Fig.[2]a). 
The swirling motions get concentrated at a height of about 500 km (which corresponds to a 
temperature minimum region, panels b and c), and then erupt (panel d). 

Because of the strong concentrations of magnetic field in the vortex core (~ 1.2 kG in 
the photosphere layer) it is interesting to consider the dynamics in terms of the evolution 
of electric currents. Figure |3] shows that the vortex tube penetration for the same event 
is accompanied by formation of an electric current sheet in the surrounding area of the 
intergranular lane (blue isosurface), which expands together with the vortex tube (yellow 
isosurface represents enstrophy) into the atmosphere and relaxes during the ejection. The 
current sheet can be a source of the Lorentz force (Sec. 5) and also, perhaps, additional 
heating. 



4. Dynamics of flow ejections 

Generally, the flow structure during the eruption phase remains twisted: the material 
around the vortex core moves up from the subsurface layers, and also towards the vortex 
from the surrounding region, and collects near the vortex edge. The plasma is moved up 
by the twisting flows into the higher atmospheric layers, and at the same time in the lower 
layers the plasma flows down though the vortex core (Fig. HJa). The magnetic field lines are 
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weakly twisted opposite to the flow direction (Fig. @]&). Because the dynamics of eruptions 
is associated with the flows surrounding a relatively narrow vortex core, we track the vortex 
core, and analyze the data in the cylindrical coordinates centered in the vortex core. For 
this analysis, we choose a typical vortex tube, and divide the vortex region in 125-km 'rings' 
(or 'zones'), where 'zone 0' includes the vortex core and its edge (Fig. [5]). 

Figure [6] shows the mean vertical velocity variations with time at different heights, from 
—300 km to 780 km, for two zones: (a) vortex core ('zone 0'), and (b) surrounding region 
('zone 1'). These diagrams show that the velocity perturbations are initiated at the depths 
of ~ 60 — 120 km below the surface, and then propagate in both directions, upward and 
downward. Such perturbations are quasiperiodic with a period of 3 — 5 min, that can be 
explained by a similar characteristic turnover time of the turbulent behavior of convection, 
where vortex tube rooted. The period is also similar to the oscillations of large-scale acoustic 
(p) modes excited in the box domain. However, these oscillations do not show correlation in 
phase with the vortex-tube oscillations, and their amplitude is essentially smaller. The mean 
vertical velocity of the perturbations increases up to 5 - 8 km/s at ~ 800 km above the solar 
surface. These upflows can be identified as small-scale flow ejections, driven by vortex tubes. 
The amplitude and quasi-period of the eruptions vary due to the vortex evolution (e.g., 
changes of their size, shape, height penetration) or/and interactions with other vortices, as 
a result of which the ejections can be magnified or suppressed. 

The comparison of the vertical velocity variations (Fig. [H]) in the vortex core (panel 
a) and the surrounding region (panel b, as indicated in Figure zones '0' and '1'), shows 
that in the vortex core the velocity perturbations have stronger amplitude. In the core, the 
upward speed of the velocity perturbations increases with height from 6 km/s in the near- 
surface layers to more than 12 km/s above 700 km. The downward perturbations propagate 
much slower, with a speed of 3 — 3.5 km/s, and their amplitude increases. In the vortex 
surrounding region (Fig. &0>), the velocity shows a similar behavior. The time shift of the 
vertical velocity variations between different zones allows us to estimate the speed of the 
flow expansion during the eruption, which is about 20 - 25 km/s in the vortex core area 
and decelerates to ~ 15 km/s at a distance of about 500 km from the vortex. For each 
individual ejection event, the estimate can vary due to the interaction of the perturbations 
with a shocks from other eruptions. 

The transformation of the velocity perturbations into shock waves additional interesting 
feature, associated with the flow ejection is (Fig. [T]c, d). This effect is identified in both the 
hydrodynamic and magnetic simulations, and in both cases associated with the vortex tube 
dynamics. In the next section we consider in more detail the process of the flow ejection and 
magnetic field effects. 
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5. Source and drivers of spontaneous flow eruptions 



The complicated dynamics of the strong swirling flows in the presence of magnetic 
field across many pressure scale-heights is an interesting interplay of hydro- and magnetic 
effects. In general, there are two type of forces that are responsible for driving the flow 
eruptions: hydrodynamic, due to pressure excess, and magnetic, caused by the Lorentz 
force. A comparison of the contributions from the hydrodynamic and magnetic effects can 
be done by estimating the RHS terms of a modified momentum equation: 



where v is the velocity vector, J is the electric current density, B is the magnetic field vector, 
p is the gas pressure, p is the density, g and c are the gravity acceleration and the speed 
of light. In this form, Eq. (1) describes the flow acceleration on the left hand side, and 
the contributions of the Lorentz force (first term) and the non-magnetic forces: the pressure 
gradient and gravity. In the initial equilibrium state the pressure gradient and gravity are 
balanced. 

Figures [7] and [8] show a comparison of the pressure gradient excess and the Lorentz 
force in the vortex core (zone '0', panels a) and the surrounding region (zone T', panels b). 
The profiles of both, non-magnetic and magnetic forces, reveal a clear connection with the 
upward and downward velocity perturbations, associated with the flow ejection, showing a 
similar decrease of the perturbation amplitude in the vortex surrounding regions, and also the 
time lag with height. Propagation of the perturbations is better visible in the Lorentz force 
(Fig. [8]) with a clear indication of acceleration in higher layers of the atmosphere, indicating 
a strong increase of magnetic field effects at the height h > 700 km. The propagation of the 
Lorentz force perturbations upward and downward along the vortex tube (Fig. [S]a) gives us an 
additional estimation for the depth of the initialization source, which varies among different 
events from the photosphere to ~ 120 km below surface. This corresponds to the non- 
magnetic case, where the primary source of the vortex eruptions is hydrodynamic. Magnetic 
field effects become important above the temperature minimum region, where upflow are 
accelerated by the Lorentz force. 

Figure M shows the contributions of the magnetic (blue curves) and nonmagnetic accel- 
eration (red curves) in four layers: 780 km and 300 km above the solar surface (panels a 
and b), the photosphere layer (panel c), and 240 km below the surface (panel d). Dashed 
curves represent the vertical velocity for the same layers, given for comparison. The results 
show that the Lorentz force is most important for the flow acceleration in the higher layers, 
where strong Lorentz-force fluctuations are correlated with the strong flow acceleration (e.g. 
at t — 9 min). The flow eruption is much weaker for the events when most contribution 
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in these layers comes from the hydrodynamic force (e.g. t = 3 min). Nevertheless, close 
to the surface and in the subsurface layers the hydrodynamic effects are significantly more 
important (Fig. M b-d) than the Lorentz force. 

As we have discussed above, the initial perturbations of the flow velocity and the hydro- 
and magnetic forces associated with the vortex eruption are generated just below the sur- 
face, where the effect of the pressure gradient force is dominant. The subsurface layers in the 
vicinity of the vortex core are characterized by strong converging flows, which compress the 
vortex tube, magnify the upward pressure gradient and accelerate the vortex-tube swirling 
motions. Figure [10] illustrates the compression process in a layer located 300 km below the 
surface (panel a), and the evolution of the vortex at the surface (panel b) with a 30-sec 
cadence. The middle snapshots in both panels correspond to the first snapshot in the atmo- 
sphere layer at the height of 625 km, presented in Figure d] and illustrate the beginning of the 
eruption. Such compression in the subsurface layers is continuous and nonstationary process, 
which is affected by convective oscillations and local turbulent magnetohydrodynamics. 



6. Conclusion 

High-resolution ground-based and space observations have revealed an intense and very 
dynamic interaction between the surface layers and the low atmosphere in quiet-Sun regions 
with relatively weak mean magnetic field. The radiative MHD simulations can reproduce 
some basic features of observed phenomena and provide an important complimentary tool 
for investigation of the underlying physical processes. Our numerical simulations with an 
initial mean 10 G vertical magnetic field show a complicated mixture of hydrodynamic and 
magnetic effects associated with spontaneous quasi-periodic (with period 3-5 min) flow 
ejections from the subsurface layers into the higher atmosphere along the magnetized vortex 
tubes, in which the magnetic field strength on the surface is typically ~ 1 kG (Fig. [6]). The 
eruptions have a complicated dynamical structure with mostly continuous swirling down- 
flows, decreasing density and heating in the vortex core, and spontaneous upflows mostly 
propagating along the vortex core periphery (Figs. [1] and |2]) and forming shock waves in the 
higher atmosphere (Figs. [T]c and d). The plasma flow in the eruptions accelerate in the higher 
(mid-chromospheric) layers from 6 to 12 — 15 km/s. Also, the perturbations, associated with 
the flow ejection, propagate into the solar interior along the vortex tube core with a speed 
about 3 — 3.5 km/s and increasing amplitude. 

The process of the flow ejection originates in a subsurface 100 km deep layer, where the 
vortex compression by converging flows increases the pressure gradient and accelerates the 
swirling flows (the evolution of a strong event is shown in Figs. [T] — |HJ) . This compressed 
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vortex tube starts penetrating into the low- atmosphere layers, and involves the surrounding 
atmospheric plasma into swirling motion. Accumulation of the swirling flows at a height of 
~ 500 km (a temperature minimum region) forms an ring-like structure, which propels the 
swirling flows into the higher layers along the vortex tube. The flows are mostly acceler- 
ated by the pressure gradient in the subsurface and near-surface layers, and by the Lorentz 
force in the higher layers (> 700 km). The described mechanism of the flow ejections in 
the vortex tubes works also in the pure hydrodynamic case, but in this case the velocity 
perturbations are much smaller and, the eruptions almost immediately fall back to the pho- 
tosphere. Nevertheless, in this case the vortex tube eruptions also generate shocks in the 
low-density atmospheric layers. 

Thus, we can conclude that the magnetic field captured in the vortex tubes plays an 
important role, magnifying the initial hydrodynamic perturbations, and accelerating them 
along the vortex tubes by the Lorentz force in the higher atmosphere, producing ubiquitous 
spontaneous flow eruptions. Our next step is to investigate the propagation of these eruptions 
into the solar corona. 
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Fig. 1. — Time-series with a cadence of 15 sec of temperature (panel a), vertical velocity (b) 
and density (c) at the height of 625 km above the solar surface; and (d) a vertical slice for 
the density, p* =log(p)) for a flow ejection initiated by a vortex tube. The location of the 
vertical slice is indicated by a short white line on the horizontal images. Time starts from 
the moment of the beginning of the vortex tube tracking for the analysis in Figures [6] - [9] 
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Fig. 2. — Different stages of the flow ejection: a) vortex tube penetration into the atmosphere 
layers, b) intensification of the swirling motions, c) concentration of the swirling motion in 
a ring-like structure, and d) flow ejection along the vortex tube. Black streamlines illustrate 
the velocity field in the vicinity of the vortex tube. Semitransparent light grey surface 
corresponds to a constant temperature of 6400 K. Yellow and blue isosurfaces correspond to 
the normalized-by-density vertical pressure gradient (— Vp/ p) for the values of 5 x 10 4 cm/s 2 
(yellow color) and —5 x 10 4 cm/s 2 (blue color). 
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Fig. 3. — Evolution of the electric current sheet during the vortex tube eruption, which drags 
it into the higher atmospheric layers (panels a-d) followed by the current relaxation during 
the flow ejection (panels e-f). Blue isosurfaces correspond to the value of electric current 
\J\ = 8 x 10 4 and yellow isosurfaces show the enstrophy distribution for Q = |curlv| 2 = 
0.35 cm -2 . Black streamlines show selected flow trajectories. 
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Fig. 4. — Topology of the flow streamlines and magnetic field lines during the flow ejection. 
The vortex is visualized by a constant temperature isosurface, T = 6000 K. The color-scale 
on this surface shows the distribution of the magnetic field strength from 50 G or less (grey) 
to 1.2 kG (red). Streamlines in panel a) illustrate the topology of flows, with helical upward 
ejection flows and downflows in the narrow vortex core. In panel b) black lines show the 
topology of the magnetic field lines. 
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Fig. 5. — Schematic representation of the analysis on zones with 125 km width around the 
vortex core. Background image illustrates the density distribution at height of 625 km above 
the solar surface for t — 7 min. For the analysis, presented in Figures [6] - EJ the vortex core 
tracked in time for the subsurface and atmosphere layers. 
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Fig. 6. — Temporal profiles of variations of the mean vertical velocity, V z in: a) vortex core, 
region '0' (range - 125 km) and b) surrounding region '1' (range 125 - 250 km) at different 
levels below the surface and in the atmosphere. Thick black curve shows the variations in 
the photosphere layer. The height difference between the curves is 60 km. 
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Fig. 7. — Temporal profiles of mean of the pressure gradient variations for: a) vortex core 
region '0' (range - 125 km) and b) surrounding region T' (range 125 - 250 km) at different 
levels below the surface and in the atmosphere. Thick black curve shows the variations in 
the photosphere layer. The height difference between the curves is 60 km. 




Fig. 8. — Temporal profiles of the mean variations of the Lorentz force normalized by density 
{Fl) z /p, for: a) vortex core region '0' (range - 125 km) and b) surrounding region '1' (range 
125 - 250 km) at different levels below the surface and in the atmosphere. Thick black curve 
shows the variations in the photosphere layer. The height difference between the curves is 
60 km. 
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Fig. 9. — Comparison of the contributions to the flow ejection of the hydro dynamic 
((— Vp/p — g)/g, red curves) and magnetic (Lorentz force, Fi,/(pg), blue curves) vertical 
accelerations (normalized by the gravity acceleration) at different altitudes: a) h = 780 km, 
b) h = 300 km, c) h = km (photosphere) and d) h = —240 km. Dashed curves correspond 
to the vertical velocity, V z , given for reference at the same layers. All curves correspond to 
the vortex core region '0'. 



-18 - 



a) 3.6 



3.4 



3.2 t 



-300 km 




3.0- 



- - - - < 



t ' ' ' 

in t t f < 1 t t 

A<\ ^ M M I 

V \ > \' ' t 1 

'Mtf t M t 
"f f I M t 

. H ' f » ' 
f f * * t I t 

tt*f * * ******* 

i,*-^ / 
•**£j Mi >-V ' 





,< /-300 km itmh 
11 / 1 

i j /,-,-km/s 
5 




2.6 2.8 



3.2 



. jr^^^jr^S * * t **^f*f f 
■J-********? / f *^* f f 
• ^s^*,*S**f f f ****** 
<*'*■*{**** $**■* — * t ** 

2.6 2.8 3 3.2 



-7 
-9 
11 



i i /0 km 
////// 



h |f ** ^ < ' - - , . , 

,* ,» * * * f **9 

****//*} 
********$, 

* * * * * * **x 
********* 
/ * * * * * / 'si 



i i i to km 7 "Tw, 

/ / / i m •* 

i { 4 i vz 
i i i jHkm/s 

t I ( i 'H r 
....... 



2.6 




2.8 3 
x, Mm 



Fig. 10. — Compression of the vortex region by surrounding converging flows. Background 
images show the vertical velocity distribution at the depth of 300 km (panel a) and in the 
photosphere layer (panel b) with cadence 30sec. Arrow show the horizontal velocity field. 
Middle snapshots in each panel correspond to the first snapshot in Figure [TJ 



